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Abstract High-latitude S 18 0 archives deriving from 
meteoric water (e.g., tree-rings and ice-cores) can provide 
valuable information on past temperature variability, but 
stationarity of temperature signals in these archives 
depends on the stability of moisture source/trajectory and 
precipitation seasonality, both of which can be affected by 
atmospheric circulation changes. A tree-ring S 18 0 record 
(AD 1780-2003) from the Mackenzie Delta is evaluated as 
a temperature proxy based on linear regression diagnostics. 
The primary source of moisture for this region is the North 
Pacific and, thus, North Pacific atmospheric circulation 
variability could potentially affect the tree-ring 8 1 ^-tem- 
perature signal. Over the instrumental period (AD 
1892-2003), tree-ring 5 18 0 explained 29 % of interannual 
variability in April-July minimum temperatures, and the 
explained variability increases substantially at lower-fre- 
quencies. A split-period calibration/verification analysis 
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found the 8 18 0-temperature relation was time-stable, 
which supported a temperature reconstruction back to AD 
1780. The stability of the 5 18 0-temperature signal indi- 
rectly implies the study region is insensitive to North 
Pacific circulation effects, since North Pacific circulation 
was not constant over the calibration period. Simulations 
from the NASA-GISS ModelE isotope-enabled general 
circulation model confirm that meteoric S 18 0 and precipi- 
tation seasonality in the study region are likely insensitive 
to North Pacific circulation effects, highlighting the paleo- 
climatic value of tree-ring and possibly other 8 lfi O records 
from this region. Our S 18 0-based temperature reconstruction 
is the first of its kind in northwestern North America, and 
one of few worldwide, and provides a long-term context for 
evaluating recent climate warming in the Mackenzie Delta 
region. 

Keywords Mackenzie Delta • White spruce tree-rings • 
Stable oxygen isotope ratios • Atmospheric circulation • 
Temperature reconstruction • NASA-GISS ModelE 

R. D. Field 

Department of Applied Physics and Applied Mathematics, 

Columbia University, New York, USA 

S. V. Kokelj 

Aboriginal Affairs and Northern Development Canada, 

Yellowknife, Canada 

T. W. D. Edwards 

Department of Earth and Environmental Sciences, 

University of Waterloo, Waterloo, Canada 


<£) Springer 


772 


T. J. Porter et al. 


1 Introduction 

Instrumental climate records for most of the planet are limited 
to the last hundred years, and the duration and spatial-density 
of these records decrease poleward (Lawrimore et al. 2011). 
Improved spatiotemporal coverage of climate data is needed 
to better characterise the high latitude climate system, assess 
its sensitivity to pre-industrial forcings, and provide initiali- 
sation and verification datasets for climate models used to 
forecast climate changes and associated impacts. In place of 
missing instrumental data, natural proxies can be used to 
fill gaps in the global climate observing network (Jones et al. 
2009; Bradley 2011). 

Tree-ring chronologies are especially useful because 
they are annually-resolved and, thus, comparable to 
instrumental climate records (Hughes 2011). The most 
climatically-sensitive tree-ring growth (e.g., ring-width or - 
density) chronologies are generally found at the margins of 
a tree’s ecological range where climate is the primary 
limiting factor (Briffa et al. 2002). Arctic treeline is one 
such environment where tree-ring growth is temperature- 
limited and has provided valuable insights on climate over 
the past millennium (Esper et al. 2002; D’Arrigo et al. 
2006; Wahl and Ammann 2007). However, in recent dec- 
ades, some white spruce trees in northwestern Canada and 
Alaska have exhibited transient temperature-growth 
responses (Jacoby and D’Arrigo 1995; D’Arrigo et al. 
2004; Wilmking et al. 2004; Pisaric et al. 2007; McGuire 
et al. 2010; Porter and Pisaric 2011), known as the 
Divergence Problem (hereafter ‘divergence’; D’Arrigo 
et al. 2008). This issue has raised concerns over the ability 
to use affected ring-width chronologies for past climate 
reconstructions (Jansen et al. 2007). 


In light of these concerns. Porter et al. (2009) explored 
the potential of tree-ring stable oxygen isotope ratios 
(§ 18 0) as an alternative proxy in the Mackenzie Delta, a 
region where divergence is widespread (Pisaric et al. 
2007). Previous studies have demonstrated potential to use 
tree-ring § 8 0 as a climate proxy (e.g., Libby and Pandolfi 
1974; Gray and Thompson 1976; Barbour et al. 2001; 
Treydte et al. 2006; Cullen and Grierson 2007), but few 
have developed tree-ring § 18 0 chronologies for long-term 
climate reconstructions (Gagen et al. 2011). Further, there 
have been few tree-ring § 8 0 studies at higher-latitudes 
(e.g., Saurer et al. 2002; Sidorova et al. 2008; Hilasvuori 
et al. 2009; Seftigen et al. 2011) where the spatiotemporal 
coverage of instrumental climate data is most limiting to 
climate change studies. 

Porter et al. (2009) developed tree-ring S 18 0 records 
(AD 1850-2003) from three white spruce trees at a site in 
the Mackenzie Delta, unofficially named ‘Timber’ (Fig. 1). 
Timber § 18 0 was shown to be most significantly associated 
with spring-summer minimum temperatures at Inuvik, the 
nearest station with a long-term climate record (AD 
1957-2003), suggesting a first-order link to temperature- 
dependent isotopic labelling of local precipitation 
(Dansgaard 1964). However, this relation has yet to be 
rigorously tested for climate reconstruction purposes. The 
main purpose of this study is to assess the robustness of the 
Timber § 18 0-temperature signal and develop an extended 
temperature reconstruction for the study region if 
warranted. 

One factor that could complicate temperature signals in 
tree-ring § 18 0 or any terrestrial 8 ls O archive that derives 
from meteoric water is variability in atmospheric circula- 
tion patterns (Edwards et al. 1996, 2008; Birks and 


Fig. 1 Map of 8 ls O study sites 
( triangles ) and climate stations 
(.solid dots ) of interest. The 
shaded relief indicates areas of 
higher (dark grey) versus lower 
(light grey) elevation 
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Edwards 2009; Liu et al. 2011; Saurer et al. 2012), which 
can alter the source region or rain-out history of atmo- 
spheric moisture, as well as precipitation seasonality, 
potentially affecting the stability of temperature signals in 
5 18 0 archives (Dansgaard 1964; Jouzel 1999; Sodemann 
et al. 2008; Sturm et al. 2010). The effects of circulation 
are most pronounced during winter and least during 
summer in the extratropics (Field 2010). Trees in the 
Mackenzie Delta do not grow during winter. However, soil 
water is a mixture of snow melt and summer rains (Mackay 
1983) and, therefore, tree-ring 5 18 0 at Timber partly 
reflects cool-season meteoric 8 18 0 and could be influenced 
by atmospheric circulation (Porter et al. 2009). 

The Mackenzie Delta is atmospherically downstream 
from the North Pacific, and the North Pacific supplies most 
of the water and energy advected to the study region and 
greater Mackenzie River Basin (Szeto et al. 2008). During 
the cold-season especially, most of this water is delivered 
by cyclonic systems from the North Pacific (Lackman and 
Gyakum 1996). As the annual Arctic sea ice melts during 
the warm season, open water in the Beaufort Sea provides 
an opportunity for ocean-atmosphere moisture flux and 
some transport inland to the study region. However, even 
during summer, the North Pacific is by far the most 
important source of moisture advected to the study region 
(Petrone and Rouse 2000). Therefore, factors that influence 
the path of North Pacific moisture to the study region 
potentially could have some effect on the Timber § ls O- 
temperature relation and could complicate the interpretation 
of the Timber S 18 0 record as a simple first-order tempera- 
ture proxy. To a large extent, the trajectory of North Pacific 
moisture is determined by the semi-permanent Aleutian 
Low (AL) over the North Pacific (Burns 1973; Wahl et al. 
1987; Dyke 2000; Serreze and Barry 2005), which varies at 
seasonal, inter-annual, decadal and even longer timescales 
(Trenberth and Hurrell 1994; Anderson et al. 2005), and is 
our main concern in this study. During times of strength- 
ened AL, atmospheric circulation becomes meridional, 
which draws warm air into northwestern North America and 
results in wetter (drier) conditions on the windward (lee- 
ward) side of the Alaska- Yukon coastal mountains (Mock 
et al. 1998; Szeto et al. 2008); a more zonal flow pattern 
occurs during periods of weakened AL. 

To determine if AL variability has influenced the Timber 
§ 18 0-temperature relation, it is important to test the relation 
over a period that experienced a wide range of circulation 
patterns. Here, we evaluate the relation using a newly com- 
piled regional temperature record, derived from six local 
datasets, that extends over the past century. Over this analysis 
period, AL has varied at interannual and decadal frequencies, 
as quantified by the North Pacific Index (NPI; Trenberth and 
Hurrell 1994). At least for winter-spring months, the range of 
NPI variability over the past century appears to be 


representative of its range over the last 400 years based on a 
well-verified NPI reconstruction by D’Arrigo et al. (2005). 
Therefore, our period of temperature-5 18 0 analysis seems 
appropriate for assessing the stability of the temperature- 
s' 8 © relation under a wide range of circulation patterns. 

However, evidence from three independent 5 I8 0 records 
in southwestern Y ukon [Fig. 1 — Mt. Logan Northwest (NW) 
Col and Prospector Russell (PR) Col ice-cores (Holdsworth 
et al. 1992; Fisher et al. 2004) and Jellybean Lake sediment- 
core (Anderson et al. 2005)] suggests a major AL regime shift 
may have occurred at ca. AD 1 840 that is unprecedented since 
at least AD 800 when a similar event is thought to have 
occurred (Fisher et al. 2004). From ca. AD 1600-1840, the 
NWCol/PRCol and Jellybean records were characterised by 
relatively high 5 18 0 values and lower 5 ls O values following 
the event. Conversely, most biological temperature proxies in 
the region indicate relatively cool conditions prior to the mid- 
nineteenth century and warmer conditions since then (Szeicz 
and MacDonald 1995; Davi et al. 2003; D’Arrigo et al. 2006; 
Viau et al. 2012). This coarse pattern of temperature vari- 
ability should have caused 8 ls O to increase in the NWCol/ 
PRCol and Jellybean records if they were dictated by a 
classical temperature effect. As such, these records were 
thought to be circulation-sensitive, and it was speculated that 
the AD 1840 shift was evidence of a shift from a zonal to 
meridional circulation regime, consistent with an intensified 
AL (Fisher et al. 2004). Simulations from an isotope-enabled 
general circulation model, however, suggest that enhanced 
meridional circulation would likely increase not decrease 
meteoric 5 18 0 (Field et al. 2010). Regardless, if a circulation 
regime shift did occur at AD 1840, it might have affected 
5 18 0 archives across interior northwestern North America 
and the stability of 5 I8 0-temperature signals in these 
archives. To explore this possibility in the Mackenzie Delta, 
we extended the Timber 8 ls O tree-ring record back to AD 
1780 to look for anomalous variability at ca. AD 1840 as 
evidence of the regime shift. 

Finally, we conducted simulations with the NASA-GISS 
ModelE isotope-enabled general circulation model 
(Schmidt et al. 2005, 2006) as an additional line of inquiry 
to better understand the potential effect of North Pacific 
atmospheric circulation on meteoric and tree-ring 5 ls O in 
the Mackenzie Delta. This exercise also provides insights 
that are pertinent to site selection for 5 ls O-based paleo- 
climate research in northwestern North America. 


2 Materials and methods 

2.1 Tree-ring 5 ls O data 

The mean climatology, site characteristics, and existing 
tree-ring 5 18 0 series for Timber site are described by Porter 
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et al. (2009). New cores from this study were sampled from 
white spruce trees in 2008 and 2009 using a Haglof increment 
borer (12 mm inner diameter). A single bark-to-bark core 
passing through or near the pith was collected from mature 
white spruce trees, providing two radii from each tree. Cores 
from three trees were selected for isotope analysis. The 
number of sample trees used in this study is low compared to 
most ring- width dendroclimatology studies. However, fewer 
independent samples are generally needed for tree-ring stable 
isotope studies to yield a robust mean site chronology (Gagen 
et al. 2011), as estimated by the expressed population signal 
(EPS; Wigley et al. 1984) or another measure of signal-to- 
noise ratio. In most tree-ring isotope studies that have 
assessed the inter-tree signal-to-noise ratio, 4-5 trees are 
often sufficient to achieve an acceptable result (Leavitt and 
Long 1984; Gagen et al. 2006), but a fixed number of samples 
cannot be prescribed for all studies due to local environmental 
conditions that can lead to more or less isotopic variability 
between trees (Leavitt 2010). In this study, only three trees 
were analysed since Porter et al. (2009) have already dem- 
onstrated that three white spruce trees yield a relatively high 
signal-to-noise ratio at Timber site. We reiterate this last point 
based on the EPS analysis discussed below. 

Sample cores were selected on the basis of tree lon- 
gevity and ring-width. Tree longevity was a criterion to 
ensure that the new tree-ring § 18 0 series would extend 
several decades prior to the AD 1840 5 18 0 anomaly 
reported by Lisher et al. (2004) in southwestern Yukon 
Territory. Ring-width was a criterion to ensure a sufficient 
sample amount for analysis. Trees at Timber have extre- 
mely narrow rings during the early-nineteenth century 
which was a technical challenge for separating individual 
rings and ensuring enough ot-cellulose after processing. To 
overcome this challenge, cores with wider tree-rings during 
the problem years were preferentially selected. This same 
criterion was used by Porter et al. (2009) due to the nar- 
rowness of tree-rings during the late-twentieth century. 

Growth patterns of the selected cores were first visually 
cross-dated, measured using a Velmex tree-ring measuring 
system, and verified against the Timber master tree-ring 
chronology using the computer program COLECHA 
(Holmes 1983). Sample rings were then separated under 
magnification using a scalpel. Corresponding rings from 
opposite radii of the same tree (i.e., representing the same 
year) were pooled, consistent with the sampling strategy 
used by Porter et al. (2009). Pooling of corresponding rings 
from the same tree is done in order to provide a more 
representative estimate of § ls O for each ring, as the § ls O 
value of a ring varies circumferentially (Leavitt 2010). This 
same sampling strategy is commonly used in traditional 
ring-width studies, but ring-width series from opposite 
radii are first measured independently and later combined 
into a mean index (Speer 2010). 


Sample rings from two of the trees span the years AD 
1790-1865, and AD 1780-1870 from the third tree. Sam- 
ples were milled to fine wood particles, approximating the 
consistency of a fine sand, and processed to '/-cellulose 
based on the method of Brendel et al. (2000) but modified 
to include an alkaline hydrolysis step to ensure complete 
removal of non-a-cellulose compounds (Gaudinski et al. 
2005; Boettger et al. 2007). In contrast, the Porter et al. 
(2009) samples were processed using a Jayme-Wise-type 
a-cellulose extraction method as described by Green 
(1963; see also Sternberg 1989). The two methods differ in 
terms of oxidants used to remove hemicelluloses. However, 
given the same wholewood, both methods produce chem- 
ically- and isotopically-identical a-cellulose when fine- 
particle samples are used (Gaudinski et al. 2005). 

The § 18 0 analysis for this study was completed at the 
G.G. Hatch Stable Isotope Laboratory, University of 
Ottawa. The system specifications and reference standards 
used are outlined by Porter and Middlestead (2012). In 
contrast, Porter et al. (2009) completed their analysis at 
the Environmental Isotope Laboratory, University of 
Waterloo. Because equipment and laboratory protocols 
used in the two studies are not identical, § 18 0 offsets 
between the two datasets can occur. A recent inter-labo- 
ratory comparison by Boettger et al. (2007) demonstrated 
this effect showing that § ls O results for the same samples 
but analysed at 9 different laboratories were systemati- 
cally offset by as much as 1 %o, and this was attributed to 
differences in sample storage and pyrolysis temperature. 
A possible analytical offset between the Porter et al. 
(2009) and new S 18 0 series was assessed by comparing 
their median chronologies over their common period of 
overlap (AD 1850-1870). An average offset of 0.7 %o was 
observed, likely due to inter-laboratory factors, and was 
corrected by adjusting the new series by +0.7 %o to bring 
them in line with the Porter et al. (2009) series. Adjusting 
the new series, however, is an arbitrary choice as it is 
difficult to say which dataset was most accurately mea- 
sured with respect to Vienna Standard Mean Ocean Water 
(VSMOW). 

To evaluate the S 18 0 precision of the samples, Quality 
Assurance (QA) standards were prepared according to 
the identical treatment method outlined by Porter and 
Middlestead (2012). Fourteen QA standards were pro- 
cessed along with 210 of the 243 samples in equal batches 
(i.e., 1 QA standard +15 samples per batch). The blanket 
5 18 0 precision estimate based on these QA standards was 
±0.3 %o, which is typical for non-processing-related 
uncertainties alone (McCarroll and Loader 2004). Given 
this result, it is likely that the samples were consistently 
processed to the same a-cellulose purity. In comparison, 
the § I8 0 precision for the Porter et al. (2009) data was 
±0.4 %o. 
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Tukey’s bi- weight robust mean is commonly used in 
dendrochronology research when combining independent 
tree-ring series into a mean chronology if sample size is 
greater than six, and the median function is used when 
sample size is less than or equal to six (Cook 1985). Both 
measures of central tendency help to guard against outlier 
data points that can strongly affect the mean chronology 
when sample size is low. Here, sample size is always less 
than or equal to six, and so the median function is used to 
calculate the master Timber 8 18 0 chronology. 

The EPS statistic (Wigley et al. 1984) was used to 
estimate the signal-to-noise ratio of the mean S ls O 
chronology: 

EPS = N x ( rbar ) 4- (1 + [N — 1) x rbar) (1) 

where rbar is the mean inter-series correlation between all 
trees over a given window length, and N is the number of 
series involved in the calculation. The EPS calculation used 
here is based on a 5 1 -year window length. Because all Timber 
5 ls O series do not cover the same period, and some series are 
not internally-continuous (e.g., null data points in the 1970s 
due to narrow rings; refer to Porter et al. 2009), the EPS 
calculation was given some flexibility to allow up to a third of 
the 5 1 years for each tree to contain null data points. In these 
cases, N is not a discrete number and was defined as the mean 
number of trees available over the 51-year period. 

2.2 Regional composite temperature records 

Regional composite minimum, mean, and maximum tem- 
perature records (AD 1892-2010) were developed and used 
to assess the long-term temperature-dependence of the 
master Timber 8 ls O chronology. The composite records 
were developed using temperature data from several local 
stations including Tuktoyaktuk (1970-2007), Inuvik 
(1957-2006), Aklavik (1926-2007), Fort McPherson 
(1892-1977), Fort Good Hope (1908-2007), and Normal 
Wells (1943-2010) (Fig. 1). Raw temperature data were 
downloaded from the Environment Canada website (http:// 
climate. weatheroffice.gc.ca). As temperatures between 
stations are biased by geographic position, the mean and 
variance of each record were adjusted to match the Inuvik 
record before calculating the monthly regional means 
(Online Resources 1-3). The regional composite records 
are discontinuous from AD 1892-1909 for most months, 
and continuous from AD 1910-2010. 

2.3 Regression model calibration and verification 

Correlation analysis was used to evaluate the relation 
between the Timber § ls O master chronology and monthly 
minimum, mean, and maximum temperatures. Monthly 
temperatures that shared a strong relation with Timber § ls O 


were averaged into an optimal seasonal index. A linear 
regression model was developed for the instrumental period 
(AD 1892-2003), with § ls O as the predictor and the optimal 
seasonal index (introduced in the Results and Discussion 
section) as the predictand. A split-period calibration/verifi- 
cation analysis (Cook and Kairiukstis 1990) was used to test 
the robustness of the temperature-8 18 0 relation. The cali- 
bration statistics were Pearson’s Product Moment Correlation 
Coefficient (r), the coefficient of determination (r 2 ), adjusted- 
r 2 . Standard Error (SE), Durbin-Watson (DW), and sign-test. 
The verification statistics were the Reduction of Error (RE) 
and Coefficient of Efficiency (CE). Upon successful valida- 
tion of the two split-period models, regression parameters for 
the full calibration period were used to reconstruct tempera- 
tures back to AD 1780 based on the Timber 8 ls O record. 

2.4 Isotope-enabled general circulation modelling 

The expected response between North Pacific atmospheric 
circulation and meteoric 8 ls O was examined using the NASA- 
GISS ModelE general circulation model, which is equipped 
with stable water isotope tracers (Schmidt et al. 2005, 2006). 
The model was configured with a 2° x 2.5° horizontal grid 
resolution and 40 vertical levels, a higher resolution than Field 
et al. (2010). Simulations were run over the 1970-2009 period 
with horizontal winds nudged toward the NCEP/NCAR 
Reanalysis (Kalnay et al. 1996) and sea surface temperatures 
and sea-ice prescribed from Rayner et al. (2003). 

Spatial correlations between NP1 (calculated as the 
mean sea-level pressure for the area 30°N-65°N, 160°E- 
140°W; Trenberth and Hurrell 1994) and the modelled 
meteoric 8 ls O were examined for fall- winter (Nov-Feb) 
and spring-summer (Apr-Jul) months. Seasonal 8 lfi O 
means, computed from monthly means, were weighted by 
monthly precipitation amount. We examined raw correla- 
tions between NPI and the meteoric 8 I8 0 field, and also 
partial correlations between NPI and a ‘residual’ meteoric 
S 18 0 field with the influence of temperature removed. Use 
of both the raw and partial correlations helped to better 
understand the extent to which any circulation influence 
could be separated from any temperature influence. Lastly, 
the relation between modelled precipitation amount and 
NPI was explored to determine the potential for changes in 
precipitation seasonality, which could influence the sea- 
sonal bias of bulk soil water 8 ls O for a given site. 

3 Results and discussion 

3.1 Timber 8 ls O chronology 

The median 8 ls O records from this study and Porter et al. 
(2009) are well correlated (r = 0.69, p < 0.001, one-tailed) 
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over the common period of overlap (AD 1850-1870), which 
indicates a common environmental forcing and supports 
combining the new and old series into a master chronology 
(Fig. 2). The master chronology spans the period AD 
1780-2003 and is defined by one to six trees depending on 
the year (overall mean is 2.9 trees/year). The average sample 
size is low compared to most traditional dendroclimatology 
studies, but EPS results indicate a high signal-to-noise ratio. 
The mean EPS for the record is 0.86 (i.e., 86 % signal, 14 % 
noise), which is above the 0.85 threshold recommended by 
Wigley et al. (1984) for dendroclimatic research. It is 
noteworthy that this relatively high EPS value was obtained 
with only three trees for most years and, to our knowledge, 
this is one of only a few tree-ring § ls O studies to demon- 
strate such a strong result with only three trees. 

In another tree-ring § l8 0 study of white spruce from the 
central Canadian sub-Arctic, Holzkamper et al. (2012) 
developed § 18 0 series from three trees and observed mean 
inter-series correlations (based on a 25-year window) 
during the period AD 1938-2003 that varied from ca. 0.3 to 
0.85, which equates to an EPS of 0.56 to 0.94. Compara- 
tively, our EPS time-series has a range of 0.74 to 0.96 over 
the last two centuries. Our EPS falls slightly below 0.85 
during most but not all years from AD 1896-1963. The 
mean EPS during this period is 0.80. The lowest EPS 
values occur in the mid-twentieth century. These lower 
values are due to the reduced sample depth during the 
1970s and a period of lower inter-series coherence during 
the 1920s and ‘30s. While much of the latter half of our 
record falls below the EPS = 0.85 threshold, these reduced 
values are only marginally less than 0.85, indicating that 
the mean chronology is still of good overall quality (always 
greater than 75 % signal, except for 1974 with an EPS of 
0.74) and warrants the use of this chronology for 


paleoclimatic interpretation. However, a higher degree 
of caution may be required when interpreting the mid- 
twentieth century portion of the chronology. 

3.2 Timber 8 18 0-temperature association 

Monthly (prior-year September to growth-year August) 
correlations between the master Timber S I8 0 chronology 
and the minimum, mean, and maximum regional composite 
temperature records were assessed over the period AD 
1892-2003 (Table 1). Consistent with the findings of 
Porter et al. (2009), we find that Timber § 18 0 is most 
closely associated with spring-summer temperatures. 
However, an exception is that we find April-July are the 
most strongly correlated months, while Porter et al. (2009) 
found March-July were most important based on their 
analysis using the shorter Inuvik temperature record (i.e., 
AD 1957-2003). Based on the longer regional temperature 
record, we find no significant correlations with March 
temperatures. 

In general, monthly correlations are strongest for 
minimum and mean temperatures. Seasonal temperature 
averages for April-May (AM), April-June (AMJ), and 
April-July (AMJJ) were correlated with Timber 5 18 0 to 
identify the optimal seasonal index. In each of these cases, 
Timber 8 18 0 was again most strongly correlated with 
minimum and mean temperatures, and we focus on these 
results. Correlations are roughly equal in strength between 
minimum and mean temperatures, and are strongest for 
AMJJ and AMJ. Split-period calibration/verification anal- 
yses were carried out to determine the seasonal index, 
minimum versus mean and AMJJ versus AMJ, that is most 
closely associated with Timber 8 18 0. The RE and sign-test 
statistics were used to identify the optimal temperature 


Fig. 2 Statistics and time-series 
for the Timber 8 ls O tree-ring 
data. Sample depth (no. trees), 
mean inter-tree correlation 
(rbar), and expressed population 
signal (EPS) are in reference to 
the individual tree 8 ls O time- 
series below; a 20-year cubic 
smoothing spline with a 50 % 
frequency cut-off (Cook and 
Peters 1981) was used to 
highlight trends in rbar and EPS. 
Individual tree 8 ls O series ( light 
grey)', median 8 ls O values of 
the Porter et al. (2009) trees 
(black) and new trees from this 
study ( medium grey)', 
interpolated 8 ls O values ( black 
dots) for the years 1972 and 
1976; master Timber 
chronology ( lower plot, black) 
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Table 1 Monthly and seasonal correlations between regional mini- 
mum, mean, and maximum temperatures and the Timber S ls O 
chronology (AD 1892-2003) 



Min 

Mean 

Max 

Year t - 1 

Sept 

0.21 

- 

- 

Oct 

0.17 

- 

- 

Nov 

0.20 

0.20 

0.20 

Dec 

0.24 

0.18 

- 

Year t 

Jan 

0.23 

0.19 

- 

Feb 

0.19 

- 

- 

Mar 

0.22 

- 

- 

Apr 

0.48 

0.46 

0.42 

May 

0.44 

0.43 

0.39 

Jun 

0.32 

0.35 

0.34 

Jul 

0.27 

0.22 

- 

Aug 

- 

- 

- 

Season 

AMJJ a 

0.55 

0.54 

0.50 

AMJ a 

0.54 

0.54 

0.51 

AM 

0.52 

0.50 

0.47 


All correlations are significant at p < 0.05 (one-tailed); correlations 
significant at p < 0.01 are bolded 


a For certain years from AD 1893-1908, temperature data (Online 
Resources 1-3) were available for April and May, and only one of 
June or July, but not both. To allow the calculation of April-July 
seasonal averages for these years, the missing June or July data points 
were estimated using the 4 closest years of data. For minimum, mean, 
and maximum temperatures, data were missing for June (1896) and 
July (1893, 1900, and 1908). Mean July (maximum July) temperature 
data were also missing for 1899 (1898 and 1899) 

index. These statistics can be thought of as tests of 
coherence, the sign-test being specific to the high-fre- 
quency domain and RE being a more general measure of 
low- and high-frequency coherence. The exhaustive details 
of these tests are not presented here. Overall, we found that 
the relation between Timber § 18 0 and AMJJ minimum 
temperatures was the strongest and most consistent 
between the split periods, and these results are presented in 
the following section. 

Last, it is worthwhile noting that April and May appear 
to contribute more to the Timber 8 I8 0 temperature signal 
than June and July given the stronger monthly correlations 
for April and May (Table 1), despite the fact that April and 
May provide only half as much precipitation to the annual 
water budget compared to June and July (28.3 vs. 54.3 cm, 
respectively) according to AD 1971-2000 climate normals 
(Porter et al. 2009). Three factors may explain this finding. 
First, April and May precipitation falls as snow, which 
tends to be swept by wind from tundra areas and frozen 
channels throughout the delta and accumulate in forested 


areas due to the snow fence effect (Rouse 1984), thereby 
likely contributing more to the bulk soil water § 18 0 upon 
melting than climate normals would suggest. Second, 
roughly half of all warm season precipitation in the 
Mackenzie River Basin is recycled from local evaporation 
and transpiration, and so only half of June-July precipita- 
tion is ‘new’ and advected from outside the region (Szeto 
et al. 2008). These two factors would tend to counter the 
apparent imbalance between April-May and June-July 
precipitation indicated by climate normals such that April- 
May would contribute more new water to bulk soil water 
than June-July. Thus, if strength of monthly temperature- 
s' 8 0 correlations is roughly proportional to the amount of 
new precipitation that falls each month, the finding that 
April and May are more important to the overall S 18 0- 
temperature signal compared to June and July would seem 
to make sense. Third, there is also the possibility that 
meteoric S 18 0 is more sensitive to temperatures during 
April-May than June-July, since the temperature-8 18 0 
relation for many parts of the world can vary seasonally 
(Fricke and O’Neil 1999). However, the lack of long-term 
meteoric S 18 0 records for this area precludes a direct 
assessment of this possibility. 

3.3 Regression model calibration and verification 

The robustness of the April-July minimum temperature-8 18 0 
relation was evaluated by a linear regression split-period 
calibration/verification analysis (Table 2). The calibration 
periods were 1893-1953 (Split 1) and 1954-2003 (Split 2), 
and the verification periods were opposite the calibration 
periods. Timber 8 ls O explains 35 and 16 % of temperature 
variability for Splits 1 and 2, respectively. The weaker 
adjusted-r 2 for Split 2 is probably due to the fact that the 
master chronology is not as well defined during the 1970s 
and ‘80s due to low sample depth; the mean chronology is 
defined by 3 trees for Split 1, and 2.2 trees for Split 2. The 
DW statistic is close to a value of 2.0 during Splits 1 and 2, 
indicating low autocorrelation in the residuals. The stan- 
dard error is lower for Split 1 (1.05 °C) compared to Split 2 
(1.36 °C), indicating more error in the Split 2 residuals, as 
would be expected given the sample depth differences. 
Results for the sign-test, a high-frequency test, are mixed. 
The Split 1 sign-test result is highly significant indicating 
that the sign of interannual 8 ls O changes closely matches 
corresponding changes in April-July minimum tempera- 
ture. However, Split 2 does not pass the sign-test. Again, 
this likely relates to the lower sample depth during Split 2. 
Despite this result, it is encouraging that Timber 8 ls O is 
associated with temperatures as strongly as it is considering 
the many potential sources of non-temperature-related 
noise, including analytical uncertainties (Werner and 
Brand 2001; Porter and Middlestead 2012), secondary 
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evaporative-enrichment effects (Roden et al. 2000; Porter 
et al. 2009), residual soil water from previous years and 
occasional flooding (Kokelj and Burn 2005), and uncer- 
tainty in the meteoric 8 8 0-temperature association 
(Dansgaard 1964; Fricke and O’Neil 1999). 

The fact that the calibration diagnostics were not exactly 
the same between Split 1 and Split 2 indicates that the 
temperature-8 18 0 relation is not perfectly stable. However, 
it is perhaps unrealistic to expect a perfectly stable relation 
given the changes in sample depth between the two periods 
and possible overprinting of the climate signal by various 
sources of analytical and environmental noise, as discussed 
above. The positive RE and CE statistics indicate that both 
of the calibration models are predictively-superior to the 
mean climatologies of the calibration and verification 
periods, respectively. This result provides some evidence 
that the temperature-8 ls O relation is relatively stable and 
suitable for reconstruction purposes. Further evidence to 
support this claim is presented below. A full-period (AD 
1893-2003) model was calibrated, explaining 29 % of 
observed temperature variability (Table 2). This level of 
temperature-8 ls O co-variability is comparable to observa- 
tions from other high-latitude tree-ring studies (e.g., 
r 2 = 0.21 to 0.47 — Sidorova et al. 2008; Csank et al. 2011; 
Seftigen et al. 2011; Holzkamper et al. 2012). A similar 
range of co-variability has been reported for other white 
spruce ring-width or -density studies in northwestern North 
America (e.g., r 2 = 0.29-0.49 — Szeicz and MacDonald 


Table 2 Calibration/verification and full-period (AD 1893-2003) 
calibration statistics for linear regression models of April-July mini- 
mum temperatures as a function of the mean Timber S ls O chronology 



Split 1 

Split 2 

Full period 

Calibration 




Period 

1893-1953 

1954-2003 

1893-2003 

N 

51 a 

50 

101 

Slope 

0.97 

0.87 

1.01 

Intercept 

-23.60 

-20.92 

-24.10 

SE 

1.05 

1.36 

1.24 

DW 

1.47 

2.19 

1.97 

r 

0.60 

0.42 

0.55 

r 2 

0.37 

0.17 

0.30 

Adjusted-r 2 

0.35 

0.16 

0.29 

Sign (+/-) 

36/14 8 

23/26 

60/40* 

Verification period 

1954-2003 

1893-1953 

n/a 

N 

50 

51 a 

n/a 

RE 

0.31 

0.46 

n/a 

CE 

0.08 

0.21 

n/a 


a The period 1893-1909 contains only 51 years of valid data points 
due to missing temperature data; sign-test results significant at 
8 p < 0.01 and * p < 0.1; all correlation results are significant at 
P < 0.01 (one-tailed) 


1995; Barber et al. 2004; Youngblut and Luckman 2008; 
Anchukaitis et al. 2012). 

The full-period calibration model is characterised by a 
DW statistic of 1.97, a standard error of 1.24 °C, and a 
sign-test of 60 agreements and 40 disagreements 
( p < 0.057) (Table 2). A more detailed analysis of resid- 
uals was used to determine if the full period model con- 
forms to standard regression model assumptions. Several 
null hypotheses were tested; (1) the residuals are not nor- 
mally distributed (Lilliefors test); (2) the residuals are not 
homoscedastic (2-sample F-test for equal variances, based 
on first and second halves of the calibration period); and (3) 
residuals are proportional to predicted values (regression 
slope test). Each null hypothesis was rejected at the 
p < 0.05 significance level, indicating that the tempera- 
ture-S 18 0 relation does not violate these critical regression 
model assumptions. 

A visual comparison of observed versus predicted tem- 
peratures (Fig. 3) provides strong evidence that the tem- 
perature- 8 18 0 relation was relatively stable during the 
instrumental period. ‘Local’ temperatures in the compari- 
son are from the regional instrumental data (Online 
Resource 1). We also used gridded temperature data from 
the Twentieth Century V2 Reanalysis (20CR) project 
(Compo et al. 2011) as additional verification. The 20CR 
data are generated from a numerical weather prediction 
model constrained by observed surface pressure, sea-sur- 
face temperatures, and sea-ice extent. For the study region 
and neighbouring areas (60-70°N, 130-140°W), the 20CR 
April-July temperatures approximate local temperatures 
reasonably well (Online Resource 4) and were used to fill 
gaps in the local temperature record from AD 1871-1909. 

Predicted and local-20CR temperatures are well mat- 
ched from AD 1871-2003 (Fig. 3). This coherence is 
especially strong at lower-frequencies (e.g., r = 0.78/0.88 
if both series are smoothed with a 5 -year/ 10-year cubic- 
smoothing spline; Cook and Peters 1981). From AD 
1871-1893, both records show a brief cooling to a low in 
AD 1893, warming to the ca. 1940s, slight cooling to the 
ca. 1960s, and a further warming trend to present. Quali- 
tatively, it is encouraging that the reconstruction captures 
the AD 1893 cooling event, but it does not fully capture the 
AD 1998 warm event coinciding with the 1997/98 El-Nino 
(Van Oldenborgh 2000) that had a major effect on global 
climate. However, given the many sources of noise that 
could influence the tree-ring 8 ls O- temperature signal, it is 
not surprising that some extreme events are not well pre- 
dicted. Overall, this comparison demonstrates a relatively 
stable temperature-8 18 0 relation. Based on the visual and 
statistical diagnostics presented, we conclude that Timber 
8 18 0 is a valid temperature proxy and use the full-period 
regression model to reconstruct April-July minimum 
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Fig. 3 Comparison between observed and reconstructed April-July 
minimum temperatures since AD 1871. Observed temperature data 
are a composite of the local and Twentieth Century Reanalysis 


(20CR) temperature data. The shaded confidence interval represents 
two standard errors of the estimate (±2.48 °C). All correlations are 
significant at p < 0.001 (one-tailed) 


temperatures since AD 1780 (Fig. 3). The reconstruction is 
discussed in the following section. 

On a final note about the relative stability of the tem- 
perature-8 18 0 relation, this finding also provides indirect 
support for the idea that Timber § I8 0 has been insensitive 
to recent North Pacific circulation variability. Over the 
twentieth century, NPI varied considerably at interannual 
and decadal timescales, especially in winter-spring months 
(Trenberth and Hurrell 1994). In spite of this variability, 
the 5 l8 0-temperature relation was found to be relatively 
stable overall (Fig. 3). Additional empirical support for this 
idea is gained by comparing April-July temperature- 8 1 8 0 
correlations between two periods defined by the upper 
(r = 0.46, 2 ct = 0.20, p < 0.05) and lower (r = 0.38, 
2cr = 0.22, p < 0.05) 50th percentiles ranked by April- 
July NPI. There is no significant difference in April-May 
temperature-8 18 0 correlations between periods of high 
versus low NPI. The implied insensitivity of Timber 8 ls O 
to North Pacific atmospheric circulation may relate to the 
spring-summer seasonal bias of Timber 8 ls O and the fact 
that extratropical meteoric S 18 0 is less sensitive to 
atmospheric circulation during the warm-season than the 
cool-season (Field 2010). Another possibility is that the 
seasonality of precipitation contributing to Timber 8 ls O is 
insensitive to North Pacific circulation effects. Both pos- 
sibilities are later evaluated using the ModelE simulations. 

Assuming that the modern range of North Pacific 
atmospheric circulation was typical of the past several 
centuries, as indicated by the D’Arrigo et al. (2005) NPI 
reconstruction, and the principle causes of variability in 
temperature and meteoric 8 ls O have been relatively 
constant, it may be reasonable to assume the Timber 8 ls O- 
temperature relation was also time-stable over the recon- 
struction period, AD 1780-2003. However, circulation 
patterns were likely very different over longer Holocene 
timescales due to changes in insolation, boreal treeline, and 
sea ice extent (Burn 1997), and these have had marked 
effects on 8 ls O archives in other regions (Edwards et al. 


1996; Hammarlund et al. 2002). As the stability of 
the Timber 8 ls O-temperature relation over these longer 
timescales is not well known, the development and 
paleoclimatic interpretation of Holocene-scale tree-ring 
S 18 0 records from this region should be guided by 
additional research on this key uncertainty. Isotope-enabled 
general circulation models forced by known changes in 
Holocene boundary conditions offer one means of 
addressing this uncertainty (LeGrande and Schmidt 2007; 
Sturm et al. 2010). 

3.4 Timber 8 18 0-temperature reconstruction 

According to the reconstruction, April-July minimum tem- 
peratures in the Mackenzie Delta region today are the 
warmest they have been in the last 230 years, which agrees 
with ring-width- and ring-density-based summer temperature 
reconstructions from Alaska, Yukon, and Northwest Terri- 
tories (Szeicz and MacDonald 1995; D’Arrigo et al. 2006; 
Anchukaitis et al. 2012), as well as multi-proxy site networks 
at the circumpolar (Overpeck et al. 1997; Kaufman et al. 
2009) and hemispheric (Mann et al. 2008) scales. The 
reconstructed mean temperature for the most recent 30-year 
period (AD 1974-2003) is —2.8 °C. By comparison, the 
coolest portion of the reconstruction occurs during the early 
nineteenth century, corresponding with the Dalton Minimum, 
a period of reduced insolation (Bard et al. 2000) and increased 
volcanism (Gao et al. 2008). The period AD 1787-1816 had 
the lowest 30-year reconstructed average of —4.6 °C. This 
cool period is well replicated by many ring-width- and ring- 
density-based temperature reconstructions across the North- 
ern Hemisphere (Briffa 2000; Luckman and Wilson 2005; 
D’Arrigo et al. 2006), and is linked to societal hardship 
(settlement abandonment, trans- Atlantic migration and war) 
in Europe (Biintgen et al. 2011). Over the 224-year 
reconstruction, April-July minimum temperatures for the 
Mackenzie Delta warmed by ca. +0.08 °C/decade or a total 
of ca. +1.8 °C. 
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3.5 The AD 1840 event 

Based on the Timber § 18 0 chronology (Fig. 2), we find no 
evidence of a major 8 18 0 shift in the Mackenzie Delta at 
ca. AD 1840. The Timber chronology exhibits a gradual 
rising trend from AD 1780 to the late-nineteenth century. 
There are no shifts from enriched to depleted values in the 
mid-nineteenth century, as were found in the NWCol/ 
PRCol (Holdsworth et al. 1992; Fisher et al. 2004) and 
Jellybean Lake (Anderson et al. 2005) records. One fluc- 
tuation in the Timber 8 18 0 record that is a departure from 
normal variability is the dip in AD 1893, but this is related 
to low temperatures recorded in the instrumental data 
(Fig. 3). Assuming Timber 5 18 0 is an archive for meteoric 
§ 18 0, these results appear to suggest that no large-scale 
North Pacific circulation regime shift occurred at AD 1840, 
or that meteoric 8 ls O and precipitation seasonality in the 
Mackenzie Delta were unaffected by the shift. 

As discussed above, the relatively stable relation 
between Timber 8 ls O and April-July minimum tempera- 
tures since AD 1871 (Fig. 3) provides indirect support for 
the idea that Timber 8 ls O is insensitive to North Pacific 
atmospheric circulation changes. We further evaluate this 
idea with the ModelE simulations, examining the influence 
of North Pacific circulation effects on meteoric S I8 0 and 
precipitation amount (Fig. 4) under the idealized condi- 
tions of a GCM. 

Results vary between November-February (NDJF) and 
April-July (AMJJ), as would be expected since North 
Pacific atmospheric circulation patterns vary widely 
between the cool- and warm-seasons (Serreze and Barry 
2005), but generalizations can be made. First we consider 
the raw NPI-8 ls O results (Fig. 4, left side). For NDJF, 
there are pronounced and inverse correlations along coastal 
Yukon near NWCol/PRCol and Jellybean Lake, indicating 
that a strengthened AL (i.e., enhanced meridional flow) 
would cause an increase in meteoric S 18 0, as was sug- 
gested by Field et al. (2010). Because the annual water 
budget for this area is cold-season biased (as per the NCEP/ 
NCAR Reanalysis; Kalnay et al. 1996), it seems likely that 
a shift in the NWCol/PRCol and Jellybean Lake 8 ls O 
records would occur given a major shift in North Pacific 
atmospheric circulation. Conversely, raw NDJF correla- 
tions are insignificant for Timber site and generally insig- 
nificant for nearby areas. However, the AMJJ season is 
more relevant to Timber 8 ls O given its spring-summer bias 
(Porter et al. 2009). The raw AMJJ NPI-8 ls O results reveal 
no significant correlation for Timber, supporting the notion 
that meteoric and tree-ring 8 ls O in the Mackenzie Delta is 
likely to be insensitive to North Pacific atmospheric cir- 
culation variability. 

Absent the effects of temperature, the partial NPI-S I8 0 
correlations (Fig. 4, center) provide important insights into 


the nature of circulation effects. For southwestern Yukon, 
the partial NDJF correlation results show weakened and 
non-significant correlations, in contrast with the raw results 
discussed above. This difference implies that much of the 
raw circulation effect for this region is mediated by tem- 
perature, consistent with a first-order Rayleigh interpreta- 
tion (Dansgaard 1964). However, some of the raw 
circulation effect along coastal Yukon cannot be accounted 
for by temperature and, thus, can be attributed to changes 
in source region and/or rain-out history. In the Mackenzie 
Delta, the partial NPI-8 ls O correlation results for NDJF 
and AMJJ show no significant correlations. Again, the 
AMJJ result is likely most relevant for identifying potential 
circulation-related complications to the Timber S I8 0-tem- 
perature signal. Overall, these results suggest meteoric and 
tree-ring 8 ls O in the Mackenzie Delta are unlikely to be 
sensitive to North Pacific atmospheric circulation 
variability. 

For Timber, NPI-precipitation correlations are generally 
weak and insignificant for NDJF and AMJJ, but with a 
greater number of significant NDJF correlations for adja- 
cent areas (Fig. 4, right; see also Online Resource 5 for a 
longer-term NPI-precipitation amount correlation analysis 
using the NCEP/NCAR [1948-2009; Kalnay et al. 1996] 
and ERA-40 [1958-2002; Uppala et al. 2005] Reanalysis 
datasets, which provides additional support for the ModelE 
results). Again, the AMJJ correlations are most important 
for Timber due to the spring-summer bias of the region’s 
water budget. Overall, the lack of any strong NPI-precip- 
itation correlations for this area suggests that precipitation 
amount and, indirectly, precipitation seasonality are 
insensitive to North Pacific circulation effects. Considering 
the NPI-precipitation and NPI-8 ls O results in tandem, it 
seems unlikely that North Pacific atmospheric circulation 
variability would strongly affect soil water or tree-ring 
S 18 0 in the Mackenzie Delta, implying that the Timber 
S 18 0 record is of considerable value as a temperature 
proxy. Given this result, the absence of the AD 1840 shift 
in the Timber 8 1!i O record should not be considered as 
evidence that a major North Pacific circulation regime 
change did not occur, as the Timber 8 18 0 record is prob- 
ably insensitive to these events. 

As for southwestern Yukon, a significant and inverse 
NDJF NPI-precipitation correlation is observed for 
NWCol/PRCol, but not Jellybean Lake, suggesting a 
stronger AL would cause increased precipitation for 
NWCol/PRCol (Fig. 4, right). However, based on the 
NCEP/NCAP Reanalysis data (Online Resources 5), NPI 
does have a significant influence on precipitation at Jelly- 
bean Lake during NDJF. No significant AMJJ NPI-pre- 
cipitation correlations are observed for NWCol/PRCol or 
Jellybean Lake. Because the annual water budget for this 
region is cold-season biased, a North Pacific circulation 
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Fig. 4 Correlation between NP1 and simulated climate fields for fall- 
winter (NDJF; top) and spring-summer (AMJJ; bottom) for the period 
1970-2009. The simulated climate fields are the raw meteoric S 18 0 
field (left), a residual meteoric S ls O field — influence of temperature 
removed (center), and the precipitation amount field (right). Corre- 
lations are significant at p < 0.05 (two-tailed; n.b., not adjusted for 


autocorrelation). Points 1, 2, and 3 are Timber, NWCol/PRCol, and 
Jellybean Lake, respectively. NPl-precipitation amount correlations 
were also assessed using the longer NCEP/NCAR (1948-2009; 
Kalnay et al. 1996) and ERA-40 (1958-2002; Uppala et al. 2005) 
reanalysis datasets (Online Resource 5) 


regime shift might be expected to induce a change in 
precipitation seasonality (i.e., ratio of cold- to warm-season 
precipitation amount). As such, the AD 1840 shift from 
enriched to depleted § 18 0 in the NWCol/PRCol and Jel- 
lybean Lake records could represent an increase in the ratio 
of cold- to warm-season precipitation amount, possibly due 
to enhanced meridional flow as was proposed by Fisher 
et al. (2004). Together, the NPI-5 18 0 and NPI-precipitation 
results indicate that the net circulation-driven response in 
§ 18 0 archives in southwestern Yukon is complex due to 
opposing effects of moisture source/trajectory and precip- 
itation seasonality, which cancel to some extent. A more 
sophisticated analysis would be required to determine 
which factor is likely to dominate in this region, but that is 
beyond the scope of this study. 

4 Concluding remarks 

Previous studies have explored the paleoclimatic potential 
of tree-ring § 18 0 for various locales around the world 
(McCarroll and Loader 2004, and references therein), but 


nearly all were limited to the modern instrumental period. 
By comparison, there are few examples of long tree-ring 
isotope chronologies developed for reconstruction pur- 
poses, a chief criticism of Gagen et al. (2011). Tree-ring 
5 18 0 has been used for quantitative climate reconstruc- 
tions, but there are few explicit temperature reconstructions 
worldwide. One example is from Ellesmere Island in the 
Canadian High Arctic where Csank et al. (2011; see also 
Ballantyne et al. 2006, 2010) used fossil tree-rings to infer 
more than 200 years of temperature variability during the 
early Pliocene. Another example is from the eastern Rocky 
Mountains of western Canada where Edwards et al. (2008) 
coupled tree-ring § 13 C and § ls O to resolve temperature and 
relative humidity signals over the past millennium. And a 
third example is from the eastern Rocky Mountains of 
Colorado where Berkelhammer and Stott (2012) recon- 
structed temperatures from tree-ring § 18 0 for the period 
AD 1485-present. Here we have demonstrated the first tree- 
ring 8 18 0-based temperature reconstruction for a high-lat- 
itude site in northwestern North America. The statistical 
and visual diagnostics of our reconstruction provide strong 
evidence that tree-ring § 18 0 is a viable climate proxy for 
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the Mackenzie Delta region and offers promise for much 
longer-term reconstructions since millennial-length tree- 
ring chronologies are known to exist in this region (Szeicz 
and MacDonald 1996; Begin et al. 2000). 

It is important to acknowledge that our reconstruction 
explains only 29 % of interannual temperature variability, 
although low-frequency trends can be more reliably pre- 
dicted. In future studies, a higher level of interannual 
variability may be explained with improved sample depth. 
Sample depth during the calibration period was especially 
low during the 1970s and ‘80s due to extremely narrow 
rings shared by most trees at Timber during this interval, 
which precluded an adequate amount of a-cellulose for 
isotope analysis since only two radii were sampled from 
each tree. The lower sample depth during this interval may 
have contributed to a lower coefficient of determination for 
the Timber 8 ls O-temperature relation. However, it is 
unclear how much more variability can be explained sim- 
ply by augmenting sample depth, as the 5 18 0-temperature 
signal is also overprinted by various sources of analytical 
and environmental noise. For future studies in the 
Mackenzie Delta and other regions where slow-growing 
trees are found, the challenge of ensuring enough ot-cellulose 
for the analysis could be easily overcome by sampling four 
or more radii from each tree and pooling corresponding 
rings into annual samples. 

The fact that the Timber § ls O-temperature signal is 
more strongly associated with spring (April-May) com- 
pared to summer (June-July) temperatures also deserves 
attention. For high-latitude boreal sites, traditional tree-ring 
indicators such as ring-width are often biased towards 
summer temperatures (Garfmkel and Brubaker 1980; 
Jacoby and D’Arrigo 1989; Briffa et al. 1998; Barber et al. 
2000; Davi et al. 2003; Wilson et al. 2007; D’Arrigo et al. 
2009; Esper et al. 2010; Andreu-Hayles et al. 2011), which 
imparts a summer bias on large-scale reconstructions that 
are based heavily on ring-width. While seasonal differ- 
ences (e.g., warm- vs. cold-season) are less important over 
multi-decadal to centennial time-scales (Jones et al. 2009), 
there is a need for proxies with a variety of seasonal biases 
in order to improve knowledge of seasonal and annual 
climate dynamics at regional- and global-scales. Improved 
knowledge of seasonal dynamics could serve several pur- 
poses including the ability to provide refined training 
datasets for climate models, and diagnosing the natural 
origins of climate anomalies. At least in the Mackenzie 
Delta, tree -ring § 18 0 responds to a broader seasonal win- 
dow of climate than tree-ring width (Begin et al. 2000; 
Pisaric et al. 2007) and, together, these proxies can 
potentially be used to tease apart differences between 
spring and summer conditions. 

The well-known relation between temperature and 
meteoric 5 18 0, especially for higher-latitude sites, presents 


an obvious opportunity to use § 18 0 archives to reconstruct 
temperatures. This relation has been most popularly 
applied to the Greenland and Antarctica ice-cores (Jouzel 
et al. 1987; Dansgaard et al. 1993; Johnsen et al. 2001; 
Augustin et al. 2004), but atmospheric circulation vari- 
ability can complicate the interpretation of § ls O archives 
as first-order proxies for temperature (Sturm et al. 2010). 
The importance of these effects for most parts of the world 
is not yet well understood, partly due to the paucity of the 
global meteoric § 18 0 records (Bowen and Revenaugh 
2003). However, as we demonstrated here, isotope-enabled 
GCMs can be used as practical tools for understanding the 
importance of circulation effects for various regions, and it 
would seem there is potential to use these models to guide 
in site selection for tree-ring § 18 0 research. For example, if 
the research objective is to develop a temperature recon- 
struction, areas where circulation effects are expected to be 
strong should be avoided. Conversely, if circulation is the 
target variable, areas that are most strongly affected by 
circulation could be sampled. 
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